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This  is  the  final  report  of  a  collaborative  research  program  by  Dr.  Jeffrey  C.  Weil  at  CIRES,  the  University 
of  Colorado,  and  Drs.  Edward  G.  Patton  and  Peter  P  Sullivan  at  the  Microscale  and  Mesoscale  Meteorology 
Division,  National  Center  for  Atmospheric  Research. 

1.  Introduction 

Dispersion  in  the  stable  boundary  layer  (SBL)  is  important  for  a  number  of  air  quality  and  other  problems,  but 
current  knowledge  of  SBL  dispersion  is  deficient  due  to  our  limited  understanding  of  the  SBL  turbulence.  This 
limitation  is  caused  by  the  weak  SBL  turbulence,  past  difficulties  in  measuring  the  turbulence,  and  the  broad 
range  of  processes  complicating  the  SBL  physics.  Over  the  past  two  decades,  large-eddy  simulation  (LES)  has 
proven  to  be  a  powerful  tool  in  simulating  turbulence  in  the  planetary  boundary  layer  (PBL).  Most  of  this  work 
has  been  aimed  at  the  convective  PBL,  but  recent  studies  have  focused  on  the  SBL  and  have  demonstrated  the 
capability  to  simulate  this  regime.  Lor  modeling  dispersion,  a  key  requirement  is  simulating  the  stochastic  nature 
of  the  turbulence,  which  LES  does  quite  well. 

Under  previous  support  from  the  Army  Research  Office  (ARO),  we  developed  Lagrangian  dispersion  models 
for  following  fluid  “particles”  in  a  turbulent  flow.  The  velocities  driving  the  particles  were  obtained  from  an  LES, 
and  the  predicted  concentration  fields  for  a  convective  PBL  were  shown  to  be  in  good  agreement  with  convection 
tank  data  ( Weil  et  al.,  2004).  Under  the  ARO  program  just  completed,  this  modeling  was  extended  to  the  SBL 
and  to  PBLs  above  a  canopy.  The  program  consisted  of  three  main  studies:  1)  dispersion  in  an  SBL  with  weak- 
to-moderate  stability;  2)  dispersion  in  a  neutral  boundary  layer  with  a  forest  canopy  at  its  base;  and  3)  LES  of  a 
dynamic  canopy  based  on  a  coupled  canopy  -  soil  model.  In  the  following,  we  discuss  the  key  results  of  these 
studies. 

2.  Dispersion  in  a  SBL 

In  a  Lagrangian  particle  dispersion  model  (LPDM),  one  tracks  passive  particles  to  obtain  the  total  dispersion  and 
mean  concentration  field.  When  using  LES  to  drive  these  models,  we  find  the  particle  velocity  from  the  LES 
“resolved”  and  “subgrid-scale”  (SGS)  velocities  at  the  particle  position.  The  SGS  velocity  is  a  random  value 
obtained  from  a  stochastic  model  using  the  SGS  turbulent  kinetic  energy. 

The  velocity  fields  used  to  drive  the  dispersion  model  were  obtained  from  the  NCAR  LES  model  ( Moeng 
and  Sullivan,  1994;  Sullivan  et  al.,  1994,  1996;  Patton  et  al.,  2005).  Two  SBL  simulations  were  run  with  the  first 
(SI)  being  a  high  resolution  version  of  the  GABLS  (Global  Energy  and  Water  Cycle  Experiment  Atmospheric 
Boundary  Layer  Study)  initiative  ( Beare  et  al,  2006).  This  LES  had  a  grid  resolution  of  ~  2  m  and  was  forced  by 
a  geostrophic  wind  speed  Ug  of  8  ms-1  and  a  surface  cooling  rate  of  0.25  Khr  The  forcings  led  to  a  surface 
friction  velocity  n*  of  ~  0.28  ms-1  and  a  stability  index  zt/L  =  1.6,  where  Zi  (=  200  m)  and  L  arc  the  SBL  height 
and  Monin-Obukhov  (MO)  length.  The  above  Zi/L  classified  the  SBL  as  weakly-stable.  The  second  case  (S2) 
was  initiated  from  the  first,  but  with  a  reduced  geostrophic  wind  (4  ms-1),  a  grid  resolution  of  ~  1  m,  a  z.i  =  62 
m,  and  a  Zi/L  =  2.2.  Lor  reference,  we  included  a  neutral  boundary  layer  (NBL,  case  Nl)  in  which  Ug  =  5  ms-1, 
iii,  =  0.25  ms-1,  and  z\  =  615  m.  Table  1  contains  a  summary  of  the  key  PBL  variables. 

These  PBLs  exhibited  similar  overall  properties  such  as  an  Ekman  spiral  with  a  mean  wind  maximum  near 
H  and  significant  wind  shear.  The  velocity  variances  had  maximum  values  near  the  surface,  and  they  decreased 
systematically  to  small  values  at  the  PBL  top  consistent  with  observations  ( Nieuwstadt ,  1984;  Lenschow  et  al., 
1987)  and  other  LESs  ( Beare  et  al,  2006).  Lor  the  SBLs,  the  Richardson  numbers  were  <  0.22  throughout  the 
boundary  layer  except  close  to  Zi-  In  addition,  the  eddy  diffusivity  for  heat  (Ku)  had  a  profile  shape  similar  to 
the  vertical  velocity  variance  g~,  profile  and  a  maximum  Ku  near  z,/ z.i  —  0.25.  The  Ku  profile  was  modeled  well 
by  Ku  =  a \lp,  where  the  lp  is  a  turbulence  length  scale  parameterized  following  Brost  and  Wyngaard  (1978). 
Overall,  the  LES  modeled  the  key  features  of  the  weakly  stable  boundary  layer  quite  well  ( Beare  et  al,  2006; 
Weil  et  al,  2006). 


Case 

R 

* 

Oi 

1 

Zi  (m) 

Zi/L  (, m ) 

SI 

0.28 

200 

1.6 

Table  1:  Flow  parameters  for  each  simulation  resulting  from 

S2 

0.13 

62 

2.2 

the  imposed  forcing  specified.  Variables  are:  «*  the  surface 

N1 

0.25 

615 

0 

friction  velocity,  Zi  the  boundary  layer  depth,  and  Zi  /L  the  bulk 

Nlc 

0.18 

612 

0 

stability  parameter  where  L  is  the  Monin-Obukhov  length. 

Calculations  of  the  mean  concentration  fields  and  dispersion  statistics  were  made  for  surface  and  elevated 
sources  in  the  above  PBLs.  The  results  were  based  on  releases  with  a  total  of  about  70,000  particles.  The  analysis 
focused  on  the  dispersion  properties  as  a  function  of  source  height  zs,  downwind  distance  x,  and  stability.  The 
following  presents:  1)  an  overview  of  the  main  dispersion  features  obtained  from  three-dimensional  (3D)  views, 
2)  vertical  and  lateral  dispersion  statistics,  and  3)  the  downwind  distribution  of  the  surface  crosswind-integrated 
concentration  (CWIC,  Cv). 

Perspective  or  3D  views  of  plumes  from  two  source  heights  in  the  SBL  reveal  key  traits  as  shown  in  Fig. 
1,  where  x  is  the  distance  along  the  PBL  averaged  wind  direction,  and  y  is  the  lateral  or  crosswind  distance. 
Near  the  source,  the  plume  has  an  anisotropic  y  -  z  cross  section  with  little  asymmetry  in  the  y  direction.  With 
increasing  distance,  there  is  an  increasing  lateral  mean  plume  displacement  due  to  the  directional  shear  or  mean 
wind  normal  to  x;  a  positive  y  deflection  occurs  at  low  elevations  with  clockwise  wind  turning  at  the  upper 
elevations.  The  directional  shear  also  leads  to  enhanced  lateral  dispersion,  lateral  asymmetry,  and  “tilting”  of  the 
plume  cross  sections  at  large  distances  (four  most  distant  planes  in  Fig.  1).  Tilting  is  a  real  feature  found  in  plume 
observations  ( Csanady ,  1973;  Smith,  1965)  especially  in  the  SBL,  and  these  results  arc  qualitatively  consistent 
with  the  observations. 

Analysis  of  the  vertical  dispersion  (cz)  and  displacement  statistics  focused  on  their  scaling.  The  vertical 
inhomogeneity  in  the  mean  wind  and  turbulence  suggested  that  their  source  height  values  were  most  relevant  for 
short-range  dispersion.  Flence,  the  az  scaled  by  is  [=  ip(zs)]  was  examined  as  a  function  of  the  dimensionless 
travel  time  or  distance  x/ (usTis),  where  £s,  us,  and  7}  v  (=  £s/aws)  arc  the  turblence  length  scale,  mean  wind  speed, 
and  Lagrangian  time  scale  at  the  source  height.  This  scaling  organized  and  collapsed  the  vertical  dispersion  results 
(az/£s)  to  approximately  a  single  curve  in  each  of  two  source  height  regimes:  zs  <  0.25 Zi  and  zs  >  0.25z,-.  The 
key  difference  between  the  regimes  was  the  different  value  (and  sign)  of  the  Kh  vertical  gradient.  Figure  2  shows 
results  for  the  lower  height  regime. 

For  reference.  Fig.  2  gives  the  statistical  theory  predictions  ( Taylor ,  1921)  applicable  to  homogeneous  tur¬ 
bulence,  and  shows  that  for  case  SI  the  simulations  arc  about  30%  lower  than  the  theory  in  the  near  field, 
x/(usTLs )  <  2.  The  difference  between  the  two  results  is  probably  due  to  the  inhomogeneity  in  the  mean  wind 
and  turbulence  in  the  simulated  PBL  and  the  grid  resolution.  For  a  finer  concentration  sampling  grid  resolution 
(factor  of  4),  the  difference  between  the  simulations  and  theory  was  only  20%.  For  case  S2,  the  dimensionless 
sampling  grid  resolution  was  coarser,  which  probably  explains  the  higher  deviation  between  the  simulations  and 
theory  in  this  case.  Overall,  the  agreement  between  the  simulated  and  theoretical  trends  is  good. 

For  the  lateral  dispersion  (ov),  the  LPDM  results  tended  towards  the  statistical  theory  predictions — <Jy  <=<  t 
and  Gy  «:  / 1/2 — at  short-  and  intermediate  times,  respectively,  and  towards  cy  f3/2  at  long  times.  The  enhanced 
spread  at  long  times  was  caused  by  the  wind  direction  shear  as  discussed  above  (for  Fig.  1)  and  led  to  plume 
tilting.  Theories  ( Saffman ,  1962;  Smith,  1965)  and  observations  have  shown  that  at  large  times  the  ay  should 
vary  as  Gy  <=<  (c)V /dz)KlJ"t3/2,  where  V  is  the  mean  crosswind  velocity  and  K-  is  the  PBL  averaged  vertical  eddy 
diffusivity.  Figure  3  shows  the  three  regimes  of  spread  for  case  SI,  where  the  long  time  (r’/2)  regime  is  clearly 
demonstrated  for  all  source  heights  and  appeal's  to  start  around  X  ~  0.3.  The  results  for  case  S2  (not  shown)  look 
similar.  In  the  NBL  (Fig.  3b),  the  spread  is  similar  although  the  /3/2  regime  appeal's  to  begin  near  X  ~  0.5,  and 
its  delay  relative  to  S 1  is  probably  due  to  the  smaller  directional  shear  in  the  NBL. 

For  the  surface  crosswind-integrated  concentrations,  a  surface  release  is  considered  first  since  it  forms  an 
upper  bound  to  the  surface  CWIC  for  elevated  sources,  and  the  results  can  be  compared  with  observations  from 
the  Prairie  Grass  experiments  ( Barad ,  1958).  Figure  4  shows  the  LPDM  predictions  for  the  SI  and  N1  cases 
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Figure  1 :  Perspective  views  of  dispersion  from  two  release  heights  in  the  SBL  with  concentration  contours  given  for  six  crosswind  (y  -  z) 
planes. 
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Figure  2:  Scaled  vertical  dispersion  as  a  function  of  the  scaled  downstream  distance  for  five  source  heights  in  the  SBL;  results  for  two 
stability  conditions  orPBLs,  SI  and  S2. 


Figure  3:  Dimensionless  crosswind  dispersion  versus  dimensionless  downwind  distance  for  five  release  heights  in  the  SBL  (SI)  and  NBL 
(Nl). 


Figure  4:  Dimensionless  cross  wind- integrated 
concentration  (CWIC)  at  the  surface  versus  the 
dimensionless  distance  for  surface  sources. 


along  with  the  Prairie  Grass  data,  where  the  scaled  CWIC  ( CyUzi/Q )  is  given  as  function  of  the  dimensionless 
distance  X  =  u,x/ ( Uzi ),  U  is  the  mean  wind  speed,  and  Q  is  the  source  strength. 

Figure  4  shows  that  the  scaled  CWIC  ( CFUzi/Q )  at  the  surface  varies  approximately  as  X  1  near  the  source 
and  as  X  1  /2  far  downstream.  The  LPDM  results  are  in  good  agreement  with  the  data  for  the  weak-to-moderate 
stability  groups:  Zi/L  =1.9  and  3.1.  For  the  more  stable  case  (Zi/L  =  5.25),  the  data  lie  above  the  curves  for 
X  >  0.07,  as  would  be  expected  due  to  the  reduced  dispersion  in  strong  stability,  which  is  not  handled  by  the 
existing  LES.  Extrapolation  of  the  LPDM  long-range  result  shows  that  the  plume  would  become  well-mixed 
(with  CyU Zi/Q=  1 )  only  very  far  downstream:  X  ~  16  or  x  ~  80  km.  This  demonstrates  the  very  slow  dispersion 
in  the  SBL. 

Papers  summarizing  these  results  were  presented  at  the  American  Meteorological  Society  (AMS)  17th  Sym¬ 
posium  on  Boundary  Layers  and  Turbulence  in  San  Diego  (May  2006)  and  at  the  Sedona  Stable  Atmospheric 
Boundary  Layer  Workshop  sponsored  by  ARO,  the  Arizona  State  University,  and  the  Defense  Threat  Reduction 
Agency  (November  2006).  A  journal  article  summarizing  the  results  is  in  preparation.  Presentations  also  were 
given  at  the  Sedona  Workshop  on  canopy  turbulence  (Patton)  and  on  subfilter  scale  motions  in  the  atmospheric 
surface  layer  (Sullivan). 

3.  Dispersion  in  a  Neutral  Boundary  Layer  with  a  Canopy 

Atmospheric  motions  interacting  with  tall  vegetation  generate  a  roughness  sublayer  (RSL)  in  which  mean  and 
turbulence  variables  depart  from  their  usual  MO  profiles  ( Finnigan ,  2000).  The  RSL  typically  occupies  the  height 
range  0  <  z/h  <  3  where  h  is  the  canopy  height.  Above  the  RSL,  the  mean  wind  and  turbulence  follow  their  MO 
similarity  forms.  Early  work  above  forests  showed  that  the  eddy  diffusivities  for  heat  and  moisture  were  about 
2  to  4  times  their  MO  values  ( Raupach ,  1979,  1981).  The  higher  diffusivities  were  caused  by  increased  scalar 
mixing  resulting  from  large  coherent  eddies  in  the  RSL  which  were  generated  by  the  inflection  point  instability 
in  the  mean  wind  at  the  canopy  top  ( Raupach  et  ah,  1996;  Finnigan,  2000). 

The  above  diffusivities  were  effectively  measures  of  dispersion  from  a  uniform  area  source.  In  this  study,  we 
focused  on  point-source  dispersion  and  investigation  of  the  height  and  downwind  extent  of  the  “canopy  effect” 
on  the  dispersion.  We  computed  dispersion  from  sources  within  and  above  the  canopy  using  the  LPDM  driven  by 
an  LES  that  explicitly  resolves  the  canopy  through  a  time-dependent  and  spatially  varying  pressure  drag  force. 
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Figure  5:  Dimensionless  crosswind-integrated  concentration  profiles  versus  dimensionless  distance  in  an  NBL  with  and  without  a  canopy. 
Upper  and  lower  panels  for  sources  within  and  above  canopy,  respectively,  where  canopy  height  (h)  is  denoted  by  dashed  line. 


For  simplicity,  the  two  LES  computations  were  conducted  for  a  neutral  boundary  layer;  one  with  a  canopy 
(case  Nlc)  and  without  (case  Nl).  The  canopy-resolving  LES  was  implemented  as  in  Patton  et  al.  (2003),  with 
the  canopy  influence  modeled  as  a  mechanical  drag  that  depended  on  leaf  area  density  and  the  square  of  the 
instantaneous  velocity.  The  canopy  was  specified  to  be  horizontally  homogeneous,  but  with  a  leaf  density  varying 
with  height,  similar  to  a  deciduous  forest  with  a  relatively  dense  over-story.  The  LES  conditions  for  the  two  NBLs 
were  essentially  the  same  except  for  inclusion  of  a  25-m  deep  canopy  in  case  Nlc. 

The  dispersion  analysis  focused  on  vertical  profiles  of  the  crosswind-integrated  concentration  (CWIC)  and 
the  particle  displacement  statistics  as  a  function  of  downstream  distance  and  source  height.  As  expected,  the 
dispersion  and  mixing  were  greater  in  the  NBL  with  a  canopy.  Figure  5  demonstrates  this  with  vertical  profiles  of 
the  dimensionless  CWIC  as  a  function  of  the  dimensionless  distance  X  for  two  source  heights:  zs/h  =  0.025  and 
3.7.  The  CWIC  profiles  are  clearly  broader  in  the  canopy  case,  and  this  was  true  for  source  heights  as  large  as 
7.5/?.  Thus,  the  canopy  effect  on  point-source  dispersion  extends  to  greater  depths  than  those  for  an  area  source. 
We  believe  this  is  due  to  the  larger  turbulence  length  and  time  scales  in  the  NBL  with  a  canopy,  a  view  that  is 
supported  by  an  earlier  dispersion  analysis  (Weil,  1990).  For  all  sources,  the  CWIC  differences  decrease  gradually 
with  X  as  a  result  of  the  rather  slow  dispersion  in  the  NBL. 

The  vertical  and  lateral  dispersion  (gv,  a.)  exhibited  some  similarities  as  well  as  differences  from  those  in  the 
SBL.  For  example,  the  vertical  dispersion  scaled  as  in  Fig.  2  collapsed  the  cz/ls  for  the  different  source  heights 
to  a  nearly  single  curve.  Also,  the  lateral  dispersion  tended  towards  the  short-  and  long-time  limits  of  statistical 
theory,  but  did  not  exhibit  the  rapid  growth  (gv  °<  t3/2)  at  large  times  due  to  directional  shear  as  found  for  the  SBL 


or  the  NBL  with  no  canopy  (Fig.  3).  This  growth-rate  reduction  was  caused  by  smaller  directional  shear  (canopy 
case),  which  was  about  50%  of  that  in  the  NBL  without  a  canopy.  The  canopy  role  in  reducing  the  directional 
shear  merits  further  investigation. 

Additionally,  the  mean  plume  height  rose  at  a  faster  rate  with  time  in  the  NBL  with  a  canopy.  This  is  attributed 
to  the  different  NBL  structural  properties  (mean  wind  and  variance  profiles)  in  the  two  cases  as  well  as  to  the 
greater  Lagrangian  time  scale  in  the  canopy  case.  A  paper  summarizing  these  results  was  presented  at  the  AMS 
18th  Symposium  on  Boundary  Layers  and  Turbulence  in  Stockholm,  Sweden  (June  2008),  and  further  analysis 
and  a  journal  article  arc  in  preparation. 

4.  LES  with  Dynamic  Canopy 

In  this  work,  we  coupled  the  LES  with  a  multi-level  canopy  model  which  solves  a  dynamic  radiation/leaf  energy 
balance  at  every  vertical  grid-level  that  will  potentially  permit  simulation  of  more  stable  boundary  layers  and 
intermittent  turbulence.  The  motivation  for  this  was  the  earlier  work  of  Van  de  Wiel  et  al.  (2002)  who  found  that 
a  canopy  layer  was  important  in  generating  intermittency  (see  Section  4.2). 

The  canopy  was  included  into  the  NOAH  (National  Center  for  Environmental  Prediction  /  Oregon  State 
University  /  Air  Force  /  Office  of  Hydrology)  land-surface  model  ( Chang  et  al,  1999),  which  describes  the 
physics  of  the  soil-water-vegetation  system.  In  addition  to  the  sensible  and  radiative  heat  fluxes,  NOAH  accounts 
for  the  moisture/vapor  flux  from  the  vegetation  and  a  soil  and  hydraulic  model  for  moisture  supply  to  the  plant 
roots.  In  its  original  form,  the  canopy  in  NOAH  exchanges  heat  and  moisture  as  a  single  “big  leaf’  and  assumes 
that  emitted  scalars  are  vented  immediately  from  the  canopy  space.  Sensible  and  latent  heat  fluxes  are  determined 
through  a  coupling  between  radiation  and  photosynthesis  models  and  the  solution  of  a  leaf  energy  balance.  In  the 
soil,  the  NOAH  model  predicts  vertical  profiles  of  temperature  and  moisture  using  a  one-dimensional  model  with 
specified  lower  boundary  conditions. 

The  key  feature  added  to  NOAH  is  a  vertically  distributed  or  “multi-level”  canopy.  NOAH  remains  as  a  1-D 
column  model  and  is  implemented  at  every  horizontal  location  in  the  LES  domain,  but  the  canopy  extends  verti¬ 
cally  into  the  flow  domain.  Therefore,  the  vertical  distribution  of  the  light  environment,  atmospheric  temperature, 
moisture  and  wind  fields  drive  the  system.  We  accomplished  this  by  merging  a  number  of  already  developed 
models  to  vertically  extend  NOAH’s  current  “big-leaf’  canopy  model. 

For  example,  the  canopy  light  environment  model  stems  from  MEGAN  ( Guenther  et  al,  1995)  which  uses 
specified  leaf  scattering  and  reflection  coefficients  combined  with  the  assumption  of  a  spherical  leaf  angle  dis¬ 
tribution  to  calculate  the  absorption  and  scattering  of  direct/diffuse  radiation  by  both  sunlit  and  shaded  leaves 
(■ Goudriaan  and  van  Laar,  1994;  Leaning  et  al. ,  1995).  Stomatal  conductance  is  calculated  using  a  photosynthesis- 
based  formulation  following  GEM  ( Niyogi  et  al,  2008)  where  a  photosynthesis  -  carbon  assimilation  -  transpira¬ 
tion  model  ( Collatz  et  al,  1991,  1992)  is  coupled  with  the  Ball-Berry  model  ( Ball  et  al,  1987;  Leaning,  1990). 

Following  Leaning  et  al  (1995),  the  ‘isothermal’  form  of  the  Penman-Monteith  combination  equation  ( Mon - 
teith,  1965;  Jones,  1976)  is  used  to  partition  net  radiation  absorbed  by  the  leaves  into  latent  heat  exchange  at  each 
level  based  upon  the  appropriate  combination  of  leaf  boundary  layer  and  stomatal  resistance  along  with  the  mois¬ 
ture  gradient  between  the  leaf  and  the  atmosphere.  Sensible  heat  exchange  between  the  leaf  and  its  surroundings 
can  be  then  determined  from  the  leaf  energy  balance  equation.  This  method  allows  for  the  heat/radiation  load  on 
the  leaf  to  determine  the  leaf’s  temperature. 

Development  of  the  multi-level  canopy  version  of  NOAH  and  its  coupling  to  the  LES  code  has  been  com¬ 
pleted.  These  models  have  undergone  two  levels  of  evaluation.  The  first  evaluation  involved  driving  the  canopy 
model  with  measured  atmospheric  conditions  taken  during  the  CHATS  campaign  ( Patton  et  al,  2008).  This 
evaluation  has  been  completed  and  showed  that  the  modeled  canopy  heat  source  distribution  matches  that  found 
from  the  CHATS  data  (Figure  6).  The  second  test  evaluates  the  coupled  model  in  its  non-canopy-resolving  mode 
to  simulate  the  effects  of  a  horizontally-varying  soil  moisture  content  on  the  turbulence  structure  of  the  convec¬ 
tive  PBL.  The  results  of  this  simulation  were  compared  with  those  of  a  previously  coupled  NOAH  -  LES  model 


Figure  6:  An  evaluation  of  the  multi-level 
canopy  model  as  driven  by  data  from  the 
Canopy  Horizontal  Array  Turbulence  Study 
(CHATS)  taken  in  a  California  walnut  orchard 
in  the  Spring  of  2007  ( Patton  et  al.,  2008).  Pre¬ 
sented  are  the  terms  in  the  vertically-integrated 
heat  budget  ( /q  ^  9'  =  SgJ  dz ),  where 

the  first  two  terms  are  taken  from  the  CHATS 
measurements  (labeled  as:  measured)  and  the 
third  (5g)  is  the  model  prediction  (labeled  as: 
modeled).  The  points  represent  hourly-averaged 
data  spanning  the  entire  last  month  of  the 
CHATS  campaign  (when  there  were  leaves  on 
the  trees).  Satisfactory  agreement  is  achieved 
using  model  parameters  (leaf  width,  length, 
scattering  and  absorbing  efficencies,  etc)  spe¬ 
cific  to  a  “broad-leaf’  forest,  and  could  be  im¬ 
proved  through  tuning  these  parameters  specifi¬ 
cally  to  walnuts. 


( Patton  et  al,  2005)  to  demonstrate  consistency  of  the  results  (not  shown). 

We  anticipate  a  third  test  (to  be  completed)  will  be  a  simulation  of  a  weakly  stable  boundary  layer  to  demon¬ 
strate  the  viability  of  the  simulated  mean  and  turbulence  profiles.  We  also  anticipate  extended  utility  from  this 
multilevel-canopy  model  through  its  inclusion  into  the  Weather  Research  and  Forecast  (WRF)  model  (NOAH  is 
already  the  primary  land-surface  model  in  WRF),  which  will  permit  incorporating  spatially  distributed  vegetation 
and  its  influence  on  turbulence  and  momentum,  heat,  moisture,  and  other  trace  gases  on  weather  and  climate 
prediction. 
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